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We investigate theoretically the dissociative ionization of a molecule using two ultrashort laser 
(pump-probe) pulses. The pump pulse prepares a dissociating nuclear wave packet on an ungerade 
surface of . Next, an UV (or XUV) probe pulse ionizes this dissociating state at large (R 
= 20 - 100 bohr) internuclear distance. We calculate the momenta distributions of protons and 
photoelectrons which show a (two-slit-like) interference structure. A general, simple interference 
formula is obtained which depends on the electron and protons momenta, as well as on the pump- 
probe delay on the pulses durations and polarizations. This interference can be interpreted as 
visualization of an electron state delocalized over the two-centres. This state is an entangled state 
of a hydrogen atom with a momentum p and a proton with an opposite momentum — p dissociating 
on the ungerade surface of . This pump-probe scheme can be used to reveal the nonlocality 
of the electron which intuitively should be localized on just one of the protons separated by the 
distance R much larger than the atomic Bohr orbit. 
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INTRODUCTION 



Recently, due to the extraordinary increase of research activities in quantum information and quantum cryptography 
there is a growing interest in various quantum intriguing phenomena originating back to the famous Einstein-Podolski- 
Rosen (EPR) paradox [l| formulated in 1935. This paradox is related to the phenomenon of quantum entanglement 
Q and to the non-local character of quantum mechanics Q or to the problem of local realism versus the completness 
of quantum mechanism. A most recent comprehensive review of various aspects of the EPR paradox can be found 
in [2|. So far nearly all experimental evidence for entanglement is related to the measurement of correlated photon 
pairs obtained in a process called " parametric down-conversion Q • In such processes a photon from a laser beam gets 
absorbed by an atom which subsequently emits two "polarization entangled" photons. The entanglement phenomenon 
should in principle also appear in various break-up processes involving slower (massive) fragments than photons, e.g. 
qq , in various disintegration processes such as photoionization and photodissociation, as suggested in So far, there 

exist very few experimental results demonstrating the entanglement in such slower processes and involving massive 
particles which are much better localized in space than massless photons Q. 

In this letter we investigate theoretically an experimental scheme based on recent advances in the ultrashort laser 
technologies which allow to shape laser pulses in femtosecond or even sub- femtosecond (attosecond) time scales [I3 . ll3l |. 
Thus these technologies allow to image the evolution in time of correlated electron-nuclear wave function. This can be 
achieved by initiating the dissociation process using an ultrashort pump pulse and allowing the dissociating fragments 
to separate and be far apart. Next, this system can be probed via a photoionization process using a probe pulse 
having a well defined phase relative to the pump pulse, and consequently in phase with the dissociating system. In 
general the photoelectron spectra will exhibit a two-centre interference, sometimes called a Fano interference since it 
was predicted by Cohen and Fano in 1967 (l4f . Their calculations showed that when a molecule is photoionized via 
absorption of one photon the photelectron spectra in diatomic molecules are modulated by an interference factor 

_ = i± sH\pe\R eq ) 
\Pe\Req 

where p e is the electron momentum and R eq is the equilibrium internuclear distance and the sign depends on the 
parity of the molecular electronic wave function; (+) for a gerade and (-) for an ungerade electronic state. We show 
that if the Fano interference is observed in dissociating diatomic molecule at large internuclear separation it may 
become an important tool for visualizing peculiarities of quantum mechanics related to entanglement and to the 
nonlocal character of the electron which (intuitively) should localize on a single heavy (1836.15 times heavier than 
the electron) centre, during the dissociation process, due to localization via Coulomb attractive force. Note that a 
quantum state describing a localized electron on a specific proton accompanying another distant proton would not lead 
to the two-centre Fano interference in the photoelectron spectrum. This interference is a result of the unique gerade 
or ungerade symmetry of the electronic molecular wave function before the photoionization takes place. Since the 
molecular Hamiltonian has also this symmetry we conclude that the state before the turn-on of the probe pulse should 
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preserve the symmetry of the dissociating state prepared by a pump pulse. In other words, quantum dissociating state 
is not a simple product of the hydrogen and proton states but is a coherent superposition of the possible " simultaneous 
presence" of the electron on both well separated protons. Thus in this specific dissociation process we face clearly the 
conflict with local realism: in any deterministic theory the simultaneous presence of the electron on two heavy well 
separated protons should not occur, and, consequently the observation of Fano interference pattern is in conflict with 
local realism in a similar way as the experimental violation Bell inequalities is 0. 
More specifically, we consider a pump-probe excitation scheme [TH, [l3[ in which a pump laser pulse prepares a 
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FIG. 1: Illustration of the proposed pump-probe experimental scheme. Two lowest electronic surfaces E 9 and E u of are 
shown. The pump pulse prepares the nuclear wave packet sliding on the upper ungerade surface E u . After the turn-off of the 
pump pulse this wave packet evolves as free system until it is photoionize at large internuclear distance R. 

dissociating nuclear wave packet on an ungerade (repulsive) surface of a molecule. Next (20-200 fs later, after 
the pump pulse is turned-off), a UV (or XUV) probe pulse ionizes this dissociating (H-atom + proton) state at 
large (R = 20 - 150 bohr) internuclear distance, as illustrated in Fig.l and also described in [ljj. We show that 
coincidence measurement of the electron and proton spectra reveal a very special, counter-intuitive nature of the 
quantum dissociation process. It can be assumed that after the turn-off of the pump pulse the motion on the ungerade 
surface of is adiabatic. Consequently, because the antisymmetric (or symmetric if disociation occurs on a gerade 
electronic state) character of the electronic wave function the quantum state of the dissociating system is very distinct 
from a state of a free proton and a free hydrogen atom [l5| which for instance occurs in proton-hydrogen scattering. 
Thus quantum mechanics predicts that even at large internuclear distance the electron can be well localized in two 
places, i.e. it can form a hydrogen atom on two well separated protons which seems counter-intuitive. Note that 
in our scheme electron localization occurs due to the Coulomb attraction from two protons. Suppose that we have 
measured the hydrogen atom at the right hand side along the laser polarization vector, as shown in Fig. 2 . Already 
at relatively small internuclear distance R=10 bohr shown in Fig. 2 the electron is very well localize on each centre. 
Thus we infer (using the charge conservation principle) that the opposite detector (at left-hand side in Fig. 2) at the 
internuclear separation larger that R=60 bohr will measure with a certitude a proton. This situation resembles a 
variant of the EPR paradox based on the disintegrating system of the two, spin one-half particles, originating from 
their initial singlet state [l6| in which by measuring a spin up by one detector we infer what was the spin projection 
measured by the opposite detector. 
Suppose now that instead of measuring directly the H-atom in the above experiment we use a probe UV pulse which 
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photoionizes the molecule and the two-centre Fano interference pattern is observed in the photoionization signal. 
Simple perturbative calculations using a plane wave approximation for the final electrons state 0, E3] predict that at 
fixed internuclear vector R the molecular photoelectron signal is modulated via interference factors sin 2 (p e • R/2) or 
cos 2 (p e • R/2). Just the fact that we observe such an interference at large \R\ would be a witness of the two following 
interesting quantum phenomena. First, in the proposed experiment the probe pulse will monitor an entangled state 
of two particles: the hydrogen atom and the proton, represented by the state \H, p > times and the H + (i,e, proton) 
flying in an opposite direction represented by a ket — p >■ The dissociating state is not a simple product of these 
basis states but because of the electron symmetry of the state dissociating on a single electron surface we must add to 
the product: \H,p> \H + , —p> the product \H, —p > \H + ,p> and integrate over the relative momentum p weighed 
with the momentum distribution of the dissociating wave packet. If a two-centre interference pattern is observed in 
the photoelectron signal this means that we observe an entangled state, since a simple product of these states cannot 
yield a two centre interference pattern. Second, any local hidden variable theory would predict the electron position 
with respect to one of the protons is well defined within the Bohr radius. Thus when the interference pattern is 
observed in the photoelectron spectrum this means that the electron is well localized simultaneously on both well 
separated protons before we photoionized this quantum system. This is a very unusual situation since the tunneling 
time from one centre to another is extremely large even at the relatively small internuclear separation R=60 bohr: 
^tunnel— 2.4 years [l^]. Moreover, if we interpret the integral of e^ 2 over half-space as a charge present around a 
specific proton one concludes that a fractional charge e/2 is well localized around one centre [18j . Thus one may 
argue that the observation of the Fano two-centre interference is a witness for the simultaneous presence of a charge 
c/2 on each centre. Clearly, this simple 3-body system with the electron and two well separated protons represents 
an interesting quantum mystery related to the formation of a hydrogen atom during the dissociation process and is 
certainly worth further experimental and theoretical investigation. 
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FIG. 2: The upper part shows schematically the first stage of the proposed experimental scheme in the case when only the laser 
pump pulse is used. This pulse dissociates the molecule. Two opposing detectors are measuring the H-atom or the H + ion 
. These measurements resemble the case of the EPR paradox in which the spin projections are measured in the disintegration 
of the two spin one-half particles from the initial single spin state. We also illustrate the definition of the Jacobi coordinates for 
the p-p-e system used in our calculations. The lower part show the two-centre Coulomb potential, the ungerade electronic wave 
function <£" ; and the corresponding probability density as function of the electron coordinate z e ; at fixed intenruclear distance 
R=10 bohr. 
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In our theoretical description of the above mentioned pump-probe scheme (dissociation followed by photoionization), 
we do not calculate the dynamics of the first step in which an ultrashort UV pump pulse photodissociates (via an 
absorption of one photon) the molecule and it thus prepares a nuclear wave packet moving on the molecular 
ungerade surface. We assume that this wave packet has the Gaussian shape right after the pump pulse is turned off, 
centres at R—Ro—12 bohr, at t=0, and it next evolves as a free system (field-free) on the £„ surface up to R=60-120 
bohr when an ultrashort UV pump laser pulse is turned-on and it photoionizes this dissociating packet also via an 
one photon process. The photoionization probability distributions of the momenta of the electron and the protons 
are calculated using first order perturbation theory in 3-D (3-D for both electron and nuclear degreees of freedom) . 
A plane wave approximation is used for the final state of protons and a "modified-plane-wave" approximation for the 
electron (see next Section for the details). Note that to the best of our knowledge, all theoretical work related to 
Fano interference has been so far done using frozen nuclei at fixed equilibrium internulear distance |i? eq |. We believe 
that our study is the first dynamical investigation of Fano interference in the photoelectron spectrum originating from 
dissociating molecules at large internuclear separations. It includes full 3-D electron-nuclear dynamics. Note, that so 
far approaches based on solving time-dependent Schrodinger equation for are based on the models with reduced 
dimensionality fl9l l20l|. 



II. PERTURB ATIVE CALCULATIONS OF DISSOCIATIVE-IONIZATION 

We use Jacobi coordinates for the electron and two protons (separated by a vector R) in which the electron position 
vector f e originates in the centre of mass of two protons, see Fig. 2. In these coordinates the Hamiltonian of H^" has 
the following form 21]: (in atomic units, h = m e = e=l) 



1 a 1 
— Af 

2m' e " 2[i 

where 



Ho = -— A r - - — + V c (f e ,R) , (2) 



2m p m e 

,H = m p /2 



2m p + m e 

are the electron - two proton reduced mass and the proton-proton reduced mass, respectively, 

V c (f e ,R) = + + -4 (3) 

\f e -R/2\ \f e + R/2\ \R\ 

is the total Coulomb interaction between protons and the electron and m e and m p are the electron and proton masses. 
The total Hamiltonian is H = Hq + V n t where 

V int = -i-^-A(t)V? e , (4) 

where 



2m p 



Vi n t describes the interaction of H^" with a UV laser probe field via its vector potential. Since we use a weak intensity 
probe pulse and the pulse frequency is larger than ionization potential of H^" we may use the perturbation theory for 
calculation of the photoionization probability of the dissociating wave packet (prepared by the pump pulse) using the 
transition amplitude: 



oo 

iHntjr „—iHot 



A fi = -i \ dt{i> f \e Mat V int e- Mot \^ in {t)) (5) 

J — oo 

We first recall the result for fixed nuclei at R. Using for the electron initial state 

ip 9 J u (7 e ,R) = ^= [pn{r c + R/2) ± i/; H (r e - R/2)] (6) 

where ipH is the hydrogen Is wave function (sign (±) is used for the electronic gerade on ungerade state), and using 
the plane wave for the final electron state: 

iP f = (2^)- 3 / 2 exp(ip e -f e ) (7) 
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we get from eq. <j5j) (just by choosing the integration coordinates local to each centre) that 

l4n 2 ~ [l ± cos(p e • R)] \A H (p e )\ 2 (8) 

where Anipe) is the atomic photoionziation amplitude given in [13]. Thus if ionization occurs from the ungerade 
surface the molecular ionization probability is modulated via a sin 2 (p e • R/2) factor. If the molecule is not initially 
aligned, and if protons momenta are not measured, we need to integrate over the direction of R which leads to the 
X factor given in eq.(pQ). To include the nuclear motion in eq.© we use the plane wave approximation for the final 
state for both the electron and for the relative motion of protons: 

ipf = (27r)~ 3 exp(ip N ■ R) exp(ip e • r e ). (9) 

We will isolate from this expression the electronic atomic photoionization amplitude Aff(p e ) in which we will use the 
exact expression for an atomic amplitude, i.e. in which exact Coulomb wave for the individual centre is included 
(but the influence of the neighboring centre is neglected which is a reasonable approximation for the case of very 
large R we are interested in). This explains the term "modified-plane-wave" approximation used in the introduction. 
Regarding the plane wave approximation for the nuclei we think that it is justified for photoionization occurring at 
very large internuclear distances. We are interested in R's as large as R > 60 bohr at which Coulomb repulsion should 
be negligible with respect to the kinetic energy of the the dissociating H^ , i.e. when 

^ » (10) 

2^ R + Po t c /n 

where Rq is the position of the the centre of the dissociating wave packet with the momentum p at t=0 when the 
pump pulse is turned off and t c is the time at which the amplitude of the probe pulse is at maximum. At this time 
the wave packet reaches the distance Rq +pot c / \i . These plane wave approximations allow us to get a simple analytic 
expression for the amplitude Afi(p e ,pN, r, AR, t c ). Using the Coulomb waves in the final state leads to much more 
complicated formulas involving some numerical integration. 

As the initial state we take the Born-Oppenheimer solution (we consider an improvement to this approximation 
in the appendix) as a product of the ungerade electronic function and the superposition of nuclear plane waves 
exp(ip • R) : 

t/jin(fe, R, Ro) = ¥ei(re,R) J d 3 p ip N {p, Rq) exp(ip ■ R) (11) 

where ipjsr(j),Ro) is the initial distribution of the momenta in the dissociating nuclear wave packet. It should be 
adjusted to the shape of the wave packet prepared by the pump pulse which we suppose is short and in the amplitude 
([5]) we assume the free evolution of the wave packet between the turn-off the pump pulse at t=0 and the turn-on 
of the probe pulse. We derive next the analytic formula for photoionization valid for any shape of <pn(p, Ro)- Its 
specific shape will be chosen for illustrating graphically our results. In order to perform analytically the time integral 
in ([5]) we need a specific shape of the vector potential A(t) of the laser field. We assume it has a Gaussian form: 



A(t) = — A exp 



(t - t ) 2 ~\ 

exp{-ifl probe t) + C.C. (12) 



2r 2 



where A , e pro f, e , r2 pro fce, T are the UV probe pulse amplitude, polarization, its central frequency and duration. Pulse 
duration r is related to the commonly used FWHM duration via relation t fwhm — 2-y/hi(2)r. We recall that thus 
defined FWHM means full width at half maximum of the laser intensity time profile, not the FWHM of the envelope 
of the laser field. t c is the time at which the probe pulse has maximum and at the same time t c is also a measure the 
time delay between the probe and the pump pulse since we have chosen t=0 as time when the pump pulse is turned 
off and the centre of the nuclear packet is at R = Rq. 

After integrating the time t and the electronic coordinate f e , in the formula ([5]) we get : 

A fi = N X A H (p e ) J d 3 p a(p) J d 3 i?exp (ip- R -ip N - i?) ( e ^ «/ 2 - e - jpVi?/2 ) (13) 

where 

a(p)=exp[if(p)t c -r 2 f(p)/2}ip N (p,R ) , (14) 
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2m' e 2(1 2\i 27r2 J /^m' e c 



A H (p e ) = Ax. 



:(Pe) = 2 3/2 (e ^ 1 e) |^ ) r ) exp \-2v arctan(p e )] j\£ , (16) 



v=— , N v = exp(7rW2)r(l + iv) . 

Pe 

Integration over the R coordinate yields two Dirac delta functions for momentum conservation. Thus we get for the 
dissociative-ionization amplitude 

A fi = N 1 {2n) 3 A H (pe) J d 3 pa(p) [(S(p-p N +p e /2) - S(p-p N -p e /2)} (17) 

Integration over the nuclear momenta p yields the final probability amplitude for dissociative-ionization 

Afi{p e ,PN> t, AR, t c ) = N 2 A H {p e ) [a(£_) - a{p+)\ (18) 

where a(p) is defined in ([T4|) . 

. ,o kA q (2it) 2 t 
2 13 ' A vn' e c 

and 

p± =pjv±p e /2 . (19) 

The shifted momentum in the last equation is related to the recoil received by each proton from the electron: the 
final relative momentum is either p + or jj_ . 

Eq.(|18j) provides a general expression for the momenta distributions of protons and of the electron valid for any 
initial distribution of momenta ippj(j),Ro) of a dissociating wave packet. However, in order to investigate in detail 
Fano two-centre interference effect we need to specify the initial momentum distribution in the wave packet prepared 
by the pump pulse. We suppose that the molecule was initially in the vibrational v=0, J=0 of the gerade bound 
state, where J is the initial angular momentum quantum number of the molecule. Thus the nuclear wave packet, after 
absorption of one photon, will be in the J=l rotational state on the ungerade E u surface of . We assume that it 
has the form: 

,_ _ . „ cosfL / AR 2 , P \ . . 

tp N (p, Ro) = C N exp I — (p-po) j exp[z(p - P)Ro\ (20) 

where cos8 p = p-e pump /\p\ and 9 p is angle between the nuclear relative momentum p and the pump pulse polarization. 
The nuclear wave packet is at R = Ro at t — 0, the probe has maximum at t = t c . This is a free Gaussian wave 
packet (in the radial variable p = \p\ ) sliding on the electronic surface £„ (its electronic wave function tp"i is given in 
© with the angular momentum J=l. The central radial momentum of the wave packet is pa and its spatial width is 
AR. In order to study the two-centre interference effect it is convenient to rewrite the probability of ionization in the 
following form (in this form the interference appears through the cross term C): 

\A fi (p* e ,p N ,T,AR,t c )\ 2 - \A H (p e )\ 2 [\a(p_\ 2 + |a(p_| 2 + C{p + , P -,t c )] (21) 

where 

(_2 A r>2 \ 

-/ 2 (p) T -^-(p-Po)J|p^/|pf, , (22) 

C{p+,p-,t c ) = 2\a(p + )\ |o(?_)|cos($(tc)) , (23) 

m c ,p e ,p N ) = (\p + \ \f4)Ro + [/(p_)| - /(?+)] t c = (\p + \ \f.\)Ro + W+{2 2 J^ 12 tc ■ (24) 
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After the use of eq. (fT9)) the phase $ becomes 

<S>(t c ,p e ,p N ) = (|p+| - |p_|) Ro+Pe ■ — t e . (25) 

A* 

This is an important result since after comparing eq.([25"j) with the phase corresponding to static result eq.® we 
conclude that by measuring the relative nuclear momentum p^ and the delay time t c we are fixing the increment of 
the intcrnuclear separation R during the time interval t c , i.e. this increment is simply: v^tc where Ujy = pjq / fj, is 
the relative velocity of protons. We can simplify more equation (1251) in the case of the nuclear momentum larger than 
the electron momentum, i.e. if the inequality p e << p^ holds we get: 



\P±\ = \/P 2 N ±Pe 'FN +Pe/ 4 — PN ± Pe • = PN ± Pe COs(6> pe ) /2; . (26) 

Thus the phase ([25]) on which the interference relies simplifies to: 

$(t c ,p N )) ~p e - R N (t c ,p N ) , (27) 

where 

R N (t c ,p N )cR -^L + ^k . (28 ) 

\PN\ A* 

We will also use later the absolute value of the vector R^ : 

R N (tc,PN) = \RN{t c ,PN)\^Ro + — ■ (29) 

M 

Clearly, we see from the last equations that, as in the case of static result (0) the interference term C(p + ,p-,t c ) 
is modulated via the term cos [RN{t c ,PN)p e cos(# ep )] , where 9 ep is the angle between the electron momentum and 
relative nuclear momentum pjy. This relation can be used for imaging the nuclear motion as suggested in |l3]: if we 
measure the ionization signal for a series of time delays t c and follow the change of a specific minimum in the spectrum, 
we can thus deduce the molecular trajectory from the relation 2rnr/p e cos(0 ep ), where n is an integer corresponding 
to a specific minimum. If furthermore, the width of the momentum distribution is sufficiently large compared to 
the electron momentum p e , i.e. p e < , we may expect that the following approximations are valid for p^ values 
close to the central value po of the momentum distributions defined via (I20|) : 



\a(p^\ ~ \a(p + \) ~\a(p N \) , (30) 



we get 



M fi | 2 ~ \A H (p e )\ 2 sm 2 [p e ■ J^LR N (t c ,\p N \)/2 



1X1 III I) 



Pn 



(31) 



\A fi \ 2 ~ cos 2 {e e )sin 2 [p e cos9 pe R N {t Cl \p N \)/2]cos 2 {e p ) (32) 

where 6 e is the angle between the electron momentum and the probe pulse polarization vector e pro £, e , and 9 p is the 
angle between the pn vector and the pump polarization e pump . Assuming that initially was at rest with the initial 
momentum of centre of mass Pcm is zero, we have the following relations between the measured final momenta 
of the two protons p\ , p-i and of the electron momentum p e and the relative nuclear momentum : 

Pi + P2 + Pe = Pcm = p N = - (pi - p 2 ) (33) 

Consequently, if the initial molecular temperature is zero it will be sufficient to measure the electron momentum and 
the momentum of one proton pi in order to determine the p/v vector on which the interference relies. Thus the vectors 
present in (1211) formulas become: 



PN = Pi + Pe/2 , p+ = pi + Pe , p- = Pi 



(34) 
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In the case of nonzero temperature of initial translational motion one should either average our formula over ther- 
mal momenta of or measure in coincidence the momenta of all three fragments resulting from the photoionization 
of dissociating in order to avoid possible washing out of the interference term. 

Summarizing, our most important result is that the Fano two-centre interference shows up in the cross term 
C(p + ,p~,t c ) in ea.([23|) via 



(fie ■ Rn^cPn) 



(35) 



where Rjsr(t c ,pisr) is given in eq. (|28p . Note, that the calculations in which R is fixed lead instead to the very similar 
interference term cos(p e • R) in eq.©. Thus the effect of nuclear motion consists in replacing R by Rn^cPn) which is 
a simple linear function of the final relative momentum of outgoing protons pjv and of the time delay t c . A convenient 
way to analyze this interference in the case of p^ fixed and parallel to the probe polarization e pro i, e (note then 9 e = 9 pe 
which simplifies significantly eq. ()32|) is to expand the angular distributions described by ()32|) in Legendre polynomials 
Pi(cos9 e ): 



\A fi 



cos 2 9 e [l — cos(p e R cos 6 e )] = cos 2 9 e 



1- J2 (2l + l)jl(PeR(tc))Pl(cOs9 e 



I— 0,1 — eve 



(36) 



\A fi \ 2 ~ /3 (t c ) +/3 2 (t c )iMcos0 e ) + Pi{t c )Pi{cos9 e ), where 



1 _ ^ sm{p e R N (t c ,p N )) 
p e R N (t c ,p N ) 



for p e R(t c ) > 1 



(37) 
(38) 
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FIG. 3: Ionization probability calculated using eqs.(21)-(24) for the case when the polarizations of the pump and probe pulses 
are parallel and the nuclear relative momentum is also parallel to both polarizations. We used: A pro f, e =60 nm, pn = Po =14.8 
a.u.,p e =0.72 a.u., AR— 3.0 bohr, 7?o = 12.0 bohr and the pump pulse duration t f hw a/ = 2.4 fs. 
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7 p e R N (t c ,p N ) 



for 7J e i?(i c ) > 1 



(39) 
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for p e R(t c ) > 1 



(40) 



We see clearly that the expansion of angular distributions in Legendre polynomials reveals the Fano interference as 
function of the pump-probe time delay t c in a very neat way. An analysis of experimental data related to the Fano 
interference using such an expansion was recently performed e.g. in [22j . More specifically, in [22j a pump-probe 
experiment was reported in which a pump pulse photodissociates a 1% molecule and the probe photoionizes the 
dissociating molecule. The ^(tdeiay) coefficient calculated from the experimental photoelectron angular distribution 
shows the modulation similar to oscillations expected from our ea.(|39p. These experimental oscillations in ^(t delay) 
do not survive for the time delays larger than few picoseconds. We suggest that this maybe related to the constant 
term present in eq. (|39|) which shows in the experiment as background or they disappear due to averaging over nuclear 
momenta which becomes more significant at larger internuclear separations. Note, that the higher j3±(t c ) coefficient 
does not contain any constant term and thus may yield a better contrast allowing the Fano interference to survive for 
larger time delays. 



III. SOME SPECIFIC EXAMPLES OF THE PROPOSED PUMP-PROBE EXPERIMENTS 



The interference expected from the theory presented in the previous section will show up most clearly when the 
proton and the electron momenta are measured in coincidence. Using our exact expressions (|2ip for probabilities as 
function of the momenta of three outgoing fragments we calculate the probabilities of ionization by the probe pulse 
for three selected geometries and plot the results in Figs. 3-6. Note that we do not use in Figs. 3-5 the approximations 
suggested in formulas ([2TH) . (150)) . The probabilities are shown as functions of the time delay t c between the pump 




FIG. 4: Same as in Fig. 3. but for the case when the polarizations of the pump and probe pulses are perpendicular and the 
nuclear relative momentum pjv is also parallel to the polarization of the as shown in the upper part of the figure. 
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and the probe pulse. We plot in Fig. 3 the angular distributions of the electron in the parallel case, i.e. all three 
vectors e pump , e pro b e , and pjv are parallel and the momenta p e and p^ are fixed at their maximal values. In Fig. 4 
the perpendicular geometry is used, i.e. we choose the case of the pump laser polarization e pump perpendicular to 
polarization of the probe pulse e pro be- In Fig. 6 again the geometry is parallel as in Fig. 5 but instead of angular 
distributions we plot there the electron spectra for the electron flying along the polarization vectors. All three 
graphs show strong interference structures as function of the time delay, as expected from the approximate factor 
sin 2 (p e ■ |^-j-i?Ar(i c ,pAr)/2^ . Note that we show in Fig. 3 and in Fig. 5 the positions of the centre of wave packet 

corresponding to certain time delays calculated using eq. (|29[) . Similar interference structures appear in the proton 
spectra displayed in Fig. 6 in which we are showing spectra as function of the single proton momentum p\ = \p± | with 
fixed electron momentum p e =0.72 a.u.. The pn vector is calculated using eq. (|34[) . 




2.18 2.20 2.22 2.24 2.26 2.28 2.30 



P e (a-u.) 



FIG. 5: Same as in Fig.3. but for the shorter pump wavelength X pro i, e — 15 nm, shorter pulse duration tfwhm=Q-'2A fs and 
smaller width of the wave packet A_R=1.0 bohr. Note that now the electron angle 8 e is fixed and equal to zero. The ionization 
probability is plotted now as function of the electron momentum p e . 



IV. CONCLUDING REMARKS 

Summarizing, we have investigated a laser pump-probe scheme in which the measurement of the momenta distribu- 
tions of the photelectron allow to determine the nuclear trajectory R(t) which simply plays a role of a slit separation in 
this double-slit like experiment in which the molecule is "illuminated from within" [13j. Thus the probe pulse prepares 
the electron source whose de Broglie wave create interference structure. The observed modulation as function of time 
delay, due to the change of the slit-separation becomes a witness of the simultaneous presence of the electron on each 
proton. In other words, when the internuclear distance R is much larger that the bohr orbit it would be natural to 
expect a localized electron on one heavy proton but such a localization would prevent the interference seen in the 
photoionization signal. 

Another unusual feature of the proposed experimental scheme in this paper is the fact that it allows to measure, 
in a sense, the sign of the spatial wave function, more specificallly the measurment we propose detects (in the case 
of an ungerade initial state) the fact that the sign of electronic wave function on one centre is minus an an another 
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remote centre is plus. Namely, at large internuclear distance the electronic wave function is given by : 

i> H (r e + R/2) ± ip H (r e - R/2) . 

Thus the probability density distribution \ip 9 J u (z)\ 2 is identical for both parities at very large R-values, see Fig. 2. 
Nevertheless, measuring the photoelectron signal allows to distinguish between the gerade and ungerade case since the 
photoionization probability at fixed internuclear distance R is proportional to cos 2 (p e • R/2) in the case of dissociation 
occurring on the gerade electronic state whereas is proportional to sin 2 (p e • R/2) in the ungerade case. This sensitivity 
of the photoelectron spectra to parity (gerade or ungerade) is very specific to this simple dissociation process and 
because of this feature this experimental scheme is very distinct form the electron two-slit diffraction. 

There exists already some experimental evidence for the existence of the Fano interferences ori gina ting from disso- 
ciating molecules at large internuclear separations: in the pump-probe experiment by Sanov et al \22\ negative iodine 
I2 ions were used in which similar to dissociating electron derealization occur when a following pump-probe is 
used. The photoionization was initialized using a 780-nm laser that dissociated the 1^ ions into I~ + neutral iodine 
atom I. After a variable time delay, a photoionizing probe removed the electron from the 1^ ion. Next the ^{ideiay) 
coefficient, defined in our eqs. (l571) .([M | was calculated from the phototelectron angular distributions. This coefficient 
oscillates as function of the time delay between the probe and pump pulses as expected from eq. (|39[) . This oscillation 
is due to the fact that as in our scheme one cannot distinguish whether the phototelectron originates from the right 
or left iodine atom separated by an internuclear distance as large as 60 bohr. 

Another method for the observation of the two-centre interference was proposed in [l5j]. This method uses the 
probe pulse which does not ionizes the dissociating molecule but is based on elastic (Thomson) photon scattering 
from the two centres in the dissociating . Thus this method, which resembles the double-slit experiment for 
photons, allows to probe the entangled state in dissociating H^" on the internuclear distance larger than the method 
discussed in our paper since the wavelength 800 nm laser is much larger than de Broglie wavelength used in our 
scheme. This method, as ours, shows the derealization of the electron on two remote centres since both the Thomson 
scattering or phoionization (in the case of our method) rely on the presence of the electron on each centre, i.e neither 



q/u 1 

Pel =-7x 




FIG. 6: Same as in Fig. 3. but now the ionization probability is plotted now as function of the proton momentum pi at fixed 
electron angle 6 e = 
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Thomson scattering n or p hotonization can occur simply on bare proton. This fact distinguishes both schemes (ours 
and that proposed in |15j ) from a simple double-slit quantum effect. 



Appendix A: Beyond the Born-Oppenheimer approximation in the initial state 

The Jacobi coordinates r e , R used in the section II are not convenient at large internuclear distances R since the 
wave functions exp(ip • R)ip^(r e ± R) used in eq.ljllj) are not the exact eigenstates of the molecular Hamiltonian even 
if one neglects the Coulomb repulsion and the attraction from the remote bare proton. These states are approximate 
eigenstates, within the Born-Oppenheimer approximation. Note, that when a proton and an hydrogen atom are very 
far apart the exact eigenstates of Ho are simply a product of two plane waves for: a free proton motion, free motion 
of the centre of mass of a hydrogen atom multiplied by the electronic wave function describing the Is electronic state 
of the hydrogen atom. Thus when the proton and the hydrogen atom are far apart and they move with relative 
momentum p it is convenient to rewrite the Hamiltonian Hq using different Jacobi coordinates from that used in 
section II. Namely, instead of using the internuclear vector R and the electron coordinate r e we now use the vectors 

-* — * — * in 

R\ = R + afi e where f\ e = f e — R/2 and a — . (Al) 

m e + m p 

Ri is the relative vector between the centre of mass of a hydrogen atom and the neighboring proton, r le is the relative 
vector between the proton and the electron. In these coordinates the Hamiltonian of H^" has the following form : (in 
atomic units, Ti = m e = e=l) 

"° " -^ A '- " ^ " W\ + Wf - S) where v > ~ + 5" (A2) 

where 

„ _ m p (m p + m e ) „ _ m p m e 

V ~ — nZ r~Z ' m e — 1 1 l A<j J 

2m p + m e m p + m e 

and m"e is simply the reduced mass of the electron in the hydrogen atom. Note that as in Jacobi coordinates used 
in section II there is no cross gradient term with coupling nuclear and electronic variables. We easily find the exact 
eigenstates of a dissociating wave packet with the relative momentum p, in these new Jacobi coordinates, in the limit 
of very large internuclear distance when the potential V2{r ei R) can be neglected. This exact asymptotic eigenstate 
has the following form: 

^in Vm (^e, R,p) = exp(ip ■ Ri)ip H (rie)) = exp(zp ■ (^R + ar le )j ^ H (f e - R/2) , (A4) 

where p is the momentum of the relative motion between the hydrogen atom and the remote proton. Note that in 
contrast to the previously used eigenstate (jlip the above state does not have the inversion symmetry with respect 
to the inversion r e — ¥ —r e . Since we expect that the initial state should have such an inversion symmetry (as being 
prepared via one photon excitation of the gerade electronic state of a H J molecule) , we construct the ungerade initial 
state in the following way 

An = 4= J d 3 ptp N (p,R ) ipZ vm (-r e ,R 7 p) - ipt° ym (r e ,R,p) . (A5) 

Clearly, for each fixed p this is an entangled state of the two particles: a free proton and a free hydrogen atom. We 
rewrite the interaction potential (j4|) in a slightly different form : 



V^t = -iA(t) ■ [ — V ?e - i— V Plp - i— V r - 2p J = -i^—A(t) ■ Vft. - i—A{t) ■ V r - 2p (A6) 
1 ™i e m p p m p P J m"c m p v 



Note that we keep here the complete interaction of the laser field A(t) with three charges whereas in the previous 
Jacobi coordinates ^ the interaction term with coupling to the system centre of mass was not included. In the last 
term in the above formula we have merged together the interaction of the electron with the proton which binds the 
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electron. The last term in eq. (|A6l will be neglected in the matrix element |5J since it describes the interaction of the 
bare proton in the case when we evaluating the term containing ipai^ie), and vice versa for the term with tpui^e)- 
Inserting this new initial state (|A5| into the eq.© we get 

A fi = N 2 [A H (p e + ap_)a(p_) - A H (p e - ap + )h(p+)} (A7) 

where 

a(p) =e X p[if(p)t c -T 2 f{p)/2}^ N (p,Ro) , (A8) 

and 

We note that using exact non-Born-Oppenheimer asymptotic (for large R) states leads to the two following modifi- 
cations in the transition amplitude as compared with (|18p. First, the atomic transition amplitude does not factorize, 
second, the function f(p) contains different reduced mass p," . The first modification has a simple interpretation: since 
the nuclei move in the opposite direction the relative electron momentum is different on each centre (proton). Since 
the atomic amplitude Ah is now different at each centre previous eqs. (|21[) - (|23[) will be modified in the following way: 

\A fi (p e ,p N ,T, AR,t c )\ 2 = \A H {pe + ap-)\ 2 \a(p-\ 2 + \A H (p e - ap+)\ 2 \a(p-\ 2 + C(p+,p-,t e ) (A10) 

where 

C(p+,p^,t c ) =2\a(p + )\ \a(p_)\\A H (p> e +ap-)\\A H (p e -ap+)\cos($(t c ,)) , (All) 



$(t c ,p e ,p N ) = (\p + \ - \p-\) R + (fHiPe + Otp-) - fH{Pe - UP+) + Pe • — *c ( Al2 ) 

A 1 

where the phase lph{p) is simply the phase of the atomic amplitude Ah defined in eg. (1161) . i.e. 

A H (p) = \A H (p)\exp[iip H {p)] . (A13) 

We conclude that the Fano interference will be similar when the non-Born-Oppenheimer correction is included. The 
only change in the time delay dependent part is a replacement of the reduced mass [i by the mass fi". Another 
change, due to shift of the argument in the Ajj{p) function modifies only the term which does not depend on the time 
delay t c . Moreover the modifications discussed in this section will be negligible for the cases studied in section III 
where the values of the electron and nuclear momenta are p e — 0.72, p^ =14.8 a.u., respectively, and a = 1/1836 is 
indeed small. Thus we do not expect that the shifts in a slowly varying function Ah (p) will modify significantly the 
predictions relative to our "dynamic" Fano interference factor (|23[) and illustrated in previous sections in Figs. 3-6. 
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